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ABSTRACT We derive the fusion hierarchy of functional equa- 
tions for critical A-D-E lattice models related to the sl{2) uni- 
tary minimal models, the parafermionic models and the super- 
symmetric models of conformal field theory and deduce the 
related TBA functional equations. The derivation uses fusion 
projectors and applies in the presence of all known integrable 
boundary conditions on the torus and cylinder. The resulting 
TBA functional equations are universal in the sense that they 
depend only on the Coxeter number of the A-D-E graph and 
are independent of the particular integrable boundary condi- 
tions. We conjecture generally that TBA functional equations 
are universal for all integrable lattice models associated with 
rational CFTs and their integrable perturbations. 



1 Introduction 

Like all good scientists, Barry McCoy has long since appreci- 
ated the power and the beauty of universality in physics and 
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its implications in mathematics. This is evident starting with 
his work on the Ising model [Q] and continues through to his 
introduction of Universal Chiral Partition Functions 0, 0]. In 
this article we follow McCoy's lead and study the universality 
of TBA functional equations. 

Ever since Baxter solved 0] the eight-vertex model, commut- 
ing transfer matrix functional equations [|6|-[TT[] have been at the 
heart of the exact solution of two-dimensional lattice models on 
a periodic lattice by Yang-Baxter methods ^ . For theories such 
as the A-D-E models considered here, these equations provide 
the key to obtaining free energies, correlation lengths and finite- 
size corrections. At criticality, the finite-size corrections are re- 
lated to the central charges and scaling dimensions of the asso- 
ciated conformal field theory (CFT). Off-criticality, these cor- 
rections yield the scaling energies of the associated (perturbed) 
integrable quantum field theory (QFT). The fundamental form 
of the functional equations involves fusion of the Boltzmann 
weights on the lattice and reflect the fusion rules of the associ- 
ated CFT. However, in order to solve for finite-size corrections 
these functional equations need to be recast in the form of a 
F-system or TBA functional equations |Tl|-p^ . Miraculously, it 
is then possible to solve for the central charges and scal- 
ing dimensions using some special tricks and dilogarithm iden- 
tities [|T| . 

More recently, it has been realized [1^, |16] that the Yang- 
Baxter methods and functional equations can be extended to 
systems in the presence of integrable boundaries on the cylin- 
der by working with double row transfer matrices. It is then 
possible to calculate surface free energies and interfacial ten- 
sions ||T^ as well as finite-size corrections and conformal parti- 
tion functions [0. The critical A-D-E models correspond, for 
different choices of regimes and/or fusion level, to unitary min- 



imal models |T3], parafermion theories [20| and superconformal 
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theories |2^. These theories include the critical Ising, tricritical 
Ising and critical 3-state Potts models. For these theories, an in- 
tegrable boundary condition on the cylinder can be constructed 
for each allowed conformal boundary condition ||2^. It is also 
possible to construct integrable seams for each conformal 
twisted boundary condition |^ on the torus. In all such cases 
it should be possible to obtain the universal conformal proper- 
ties in the continuum scaling limit by solving suitable functional 
equations. 

In this paper we derive general fusion and TBA functional 
equations for the critical A-D-E lattice models. Although the 
fusion hierarchy of functional equations is not universal, we show 
in this paper that the F-system or TBA functional equations for 
the A-D-E models are universal in the sense that they depend 
on the A-D-E graph only through its Coxeter number, and more 
importantly, they are independent of the choice of integrable 
boundary conditions. 

The universality of the TBA equations has important conse- 
quences. It asserts that the functional equations are the same 
for all twisted boundaries on the torus and open boundaries 
on the cylinder. Therefore the same functional equations must 
be solved in all cases! So no new miracles, beyond the periodic 
case, are required to solve these equations in the presence of 
conformal boundaries. Instead, the different solutions required 
among the infinite number of possible solutions to the TBA 
functional equations are selected by appropriate analyticity re- 
quirements. These analyticity properties allow for the deriva- 
tion of non-linear integral equations (NLIE) that can be solved 
for the complete spectra of the transfer matrices and the uni- 
versal conformal data encoded in the finite-size corrections. Of 
course the analyticity properties are not universal. However, one 
strength of the lattice approach is that the analyticity deter- 
mined by the structure of zeros and poles of the eigenvalues of 
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the transfer matrices can be probed directly by numerical cal- 
culations on finite-size lattices. In this way it is possible to build 
up case by case a complete picture of the required analyticity 
properties. 

The layout of the paper is as follows. We first recall some 
results about fused A-D-E models in Sections lA-'LS. In Sec- 
tion ^, we define the transfer matrix for the different boundary 
conditions, on the torus and on the cylinder, with and without 
seams. In Section ^, we state the main result of the paper, that 
is the TBA equation, the boundary specific functional equa- 
tions and their universal form. In Section ^, we derive the TBA 
and related functional equations. We first study the general idea 
which is based on local properties in and we then apply it 



to the torus in |4.2| and the cylinder in f4.3| . We conclude with a 
discussion in Section ^. 

The methods developed in this paper should extend to the 
general si{2) coset models. However, we focus our attention on 
the unitary minimal, superconformal and parafermionic series, 
corresponding to fusion levels p = 1,2 and negative regime re- 
spectively, because only in these cases is our knowledge of the 
integrable and conformal boundary conditions complete. 

Our results can also be extended, using the methods of []TB[, to 
the A and D lattice models off-criticality yielding precisely the 
same TBA functional equations. In these cases the lattice mod- 
els admit elliptic solutions to the Yang-Baxter equation where 
the elliptic nome plays the role of the deviation from critical 
temperature. Although the TBA functional equations are tradi- 
tionally associated with the thermodynamics of quantum spin 
chains at finite temperature T they are derived here from a two- 
dimensional lattice approach. It is of course well known that the 
relevant quantum spin chains can be obtained as a logarithmic 
derivative of the transfer matrices with respect to the spectral 
parameter. In fact, the two approaches based on spin chains and 
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two-dimensional lattice models are entirely equivalent p, |25[| . 
Most importantly, the TBA equations originally conjectured by 
Zamolodchikov []T^ can be derived within the lattice ap- 
proach yielding the precise relation between the temperature 
T and the elliptic nome of the lattice model. The critical case 
that we focus on here just corresponds to T = and vanishing 
elliptic nome. 

It is worthwhile to mention possible applications of the TBA 
functional equations. Given prescribed conformal boundary con- 
ditions, the TBA functional equations can of course be used to 
derive non-linear integral equations (NLIE) which in principle 
can be solved for the known conformal partition functions. Of 
more interest, however, is to include in the NLIE the effect of 
a boundary field C, to perturb away from the conformal bound- 
ary. This induces an integrable boundary flow and allows for the 
study of integrable boundary flows between distinct conformal 
boundary conditions [Q. In the off-critical case for the A and 
D models, a similar analysis [^] allows for the study in the 



presence of boundaries of thermal renormalization group flows 
connecting different coset models. 

1.1 Face Weights 

A lattice model in the A-D-E series is associated with a graph 
G, of A, D oi E type. The spins are nodes of the graph G and 
neighbouring sites on the lattice must be neighbouring nodes of 
the graph. The probability distribution of spins is defined by the 
critical (unfused) Boltzmann weight of each face (or plaquette) 
of spins, depending on a spectral parameter u: 

d c I 

Wi'^ I v] = =s{\-u)5^, + s{u)J'l^5,, 

\a b J [\ I y 

' (1.1) 

where g is the Coxeter number of G, A = ^, s{u) = and ipa 

is the entry, associated with the node a, of the Perron- Frobenius 
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eigenvector of the adjacency matrix G. 

These Boltzmann weights are represented by a local face op- 
erator Xj{u) in the Temperley-Lieb algebra T{N, A) |T^: 




s{X — u)l + s{u)ej 



[1.2) 



i-i j i+i 



where = Xj{X) is a Temperley-Lieb generator and j = 1, . . . , N 
labels the position in the lattice. 



1.2 Fusion Projector 

In turn, this model gives rise to a hierarchy of fused models 
whose Boltzmann weights we are going to describe. 

We first define recursively the fusion operators PJ, for r G 
(1, g) as follows: 



Pi 



pr 

j 



(r-2)A)P;;^ r>3, 



;i.3) 



where Sk = s{kX) and j is suitably restricted [|T9|. Thus, PJ can 
be expressed as a function of Cj, Cj+i, . . . , Cj+^r-s)- In particular. 




>-'2 




;i-4) 



We shall represent the fusion operators diagrammatically as 



pr 

j 



:i.5) 



j+r-3 j-+r~2 



It is easy to show that this operator is in fact a projector. 
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Moreover, 



pr' pr _ pr pr' _ pr 
^3' ^3 - ^3 ^3' - ^3 



for 0<j' — j<r — r'. 



(1.6) 



This implies that any operator expressible as a product of local 
face operators and falling within the boundaries of a projector, 
acts as a scalar on it: 

Xj,{u)P^ = PJXj,{u) = s{X - u)PJ for < / - i < r - 2. 

(1.7) 



s(A — u) 



A particularly important case is for u = +A: its local face oper- 
ator is a projector orthogonal to all the PJ. 

As Pj is clearly translationally covariant (in its domain of 
definition) we can decompose it onto the spaces of paths with 
given end points: P^{a, b) is the fusion projector acting on paths 
from a to 6 in r— 1 steps. Its rank is given by the fused adjacency 
matrix entries: 

Rank (P^(a,6))=F:, (1.8) 

also called basic intertwiners and recursively defined by the si2 
fusion rules: 



= 1, = G, F'" = F'^-^F^ - F^-^ for r = 3, . . . , c/ 

This equation can be recast as 

(F")2 = (/ + ^'-i)(/ + F''+^) 



(1.9) 
1.10) 
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where 

pr ^ pr-l pr+l _ (l-H) 

The +1 eigenvectors of P''(a, 6) are thus indexed by an integer 
7 G (1, F^^) refered to as the bond variable. We denote them by 
C/^(a, b) and call them fusion vectors. 



1.3 Fused face operators 

These projectors allow us to define the (p, g)-fused face operator 
defined as the product of q rows of p local face operators with a 
shift of the spectral parameter by ±A from one face to the next: 




(1.12) 



The position of the projectors and spectral parameters can be 
altered by pushing-through: 




These properties imply several others, namely the Transposition 
Symmetry, 

Xf{uY = Xf{u + {q-p)\l (1.14) 
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(1.15) 



(1.16) 



where s\'^{u) — 11^=0 11 fc=o •^('^ + ^"^^ the Abelian 

Property, 



(« + {p-l)X)Xf{v + (g-l)A) = X^'iv + {p-l)X)Xfiu + {q- 

(1.17) 



These operators, contracted against the fusion vectors, yield 
the {p, g)-fused Boltzmann weights. They depend not only on 
the spins on the four corners but also on bond variables on the 
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edges: 



( d 


7 


c 














a 


b 





up — 




U „ (0,6) 



(1.18) 

where the normahzation function Sq^^(u) eUminates some scalar 
factors common to all the spin configurations which appear in 
the process of fusion. In the case, the bond variables are 
trivial, that is, a, /5, 7, 5 = 1. 

The fused Boltzmann weights satisfy a Diagonal Reflection 



/ d 


7 


c 














a 


b 





d 6 a 



{u) 



W^Pl 7 



a 



c P b 



u + {q-p)X 
(1.19) 



and Crossing Symmetry: 



d 'J c 
WP'^i 6 13 
a a b 



U 



u) 



W 



IP 



a 


S 


d 




a 




7 


A 


b 




c 





X-i 
(1.20) 



2 Transfer matrix 

Given the hierarchy of fused Boltzmann weights, we build trans- 
fer matrices for different fusion levels and boundary conditions: 
on the torus, and on the cylinder, with or without seams. 



2. 1 Seams 

Simple scams arc modified faces. They come in three different 
types, r, s and C-type. A label (r, s, Q G x Ag-i x F, where 
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r is the symmetry algebra of the graph G, encodes a triple seam 
involving three modified faces. The symmetry ( is taken as the 
identity when omitted. 

We first define the r-type seam for the (p, g)-fused 

model. It is a usual (r — 1, g)-fused face (it does not depend 
on the horizontal fusion level p) with a boundary field C, acting 
as a shift in the spectral parameter, and another choice for the 
removal of the common scalar factors: 



/d 7 


c 








1 — — 1 




P 






V 

\a a 


b 











^r— 1 q—1 



r-2qi 



An s-type seam is the normalized braid limit of an r-type 
seam, it does not depend on any spectral parameter: 

d ^ c 




(1,^)1 



.■(9+i)(«-i)g - 



lim 

>joo 



Mi 



( d ^ c 




s p 




\ a a b 




(2.1) 




an r-type 




/ d •y c 




s p 




y a a b 





(2.2) 

The automorphisms C G F of the adjacency matrix, satisfying 
Ga,b = G({a),cib), leave the face weights invariant 

fCid) 7 C(c) 



1 d 


7 


c 








P 




V 

\ a 


a 


b 





d cf c 



u /i 



u /( 



C{a) 



Vc(c 



a Cib) 
(2.3) 



u 



and act through the special seam |28 
d c \ d ad) 




Sb<:{a)Scad) ^< f 



C(a) 



1, 

0, 



F^d' ^ 0, a = /? 



Notice that the (r, s, () = (1, 1, 1) seam, where ( 
the identity automorphism, is the empty seam 

d c 

\ P \ = Sab Scd 6afS F^t^ 



g+1 

b = ({a), c = C{d), 
otherwise. 

(2.4) 
1 denotes 



(2.5) 
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The push-through property is also trivially verified for a C-type 
seam. 

The label s appearing in a (1, s)-seam is an integer in Ag_i. 
In ||2^, we define an (l,a)-seam with a E G which reduces to 



the definition given here for G oi A type but which extends it 
for the -Deven, Eq and graphs. 



2.2 Torus transfer matrix 

The transfer matrix for the (p, g)-fused model with an (r, s, Q- 
seam, on the N faces torus on the square lattice is given, in 
the basis of the cychc paths in N steps plus the seam, with 
bond variables between adjacent spins, by the product of the 
corresponding Boltzmann weights: The entries of the transfer 
matrix with an (r, s, C) seam are given by 



b2 



-ftv 1 — PN+l 1 — PN+2 



0-N 



'■(",0 



111 



'N \ l 



ajv+2 



N 



7 





where the sum is over all possible vertical bond variables. The 
usual periodic boundary condition is obtained for (r, s, Q = 
(1, 1, 1). The definition can be generalised to accommodate an 
arbitrary number of seams. Because the seam faces are modi- 
fied bulk faces, they satisfy the GYBE, so they can be moved 
around freely with respect to the bulk faces, the spectrum of 
the corresponding transfer matrices remains unchanged. How- 
ever, in the D2k cases, when there are several seams, their order 
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can not be exchanged because the fusion algebra of defect lines 
is non commutative PD[| . 



2.3 Boundary weights 



The boundary weights are labelled by (r, a) with r G Ag_2 a 
fusion level and a e G a node of the graph. 

In the Al case, all (r, s) boundary weights are obtained from 



the action of an (r, s)-seam on the vacuum boundary weight 
and we construct in an (r, a)-seam with a G G so that it is 
also the case for the -Dgven, Eq and Eg graphs. Nevertheless, in 
all cases, the (1, a) boundary weights, for two g-adjacent nodes 
of G, c and a (i.e. -F^jt^ 7^ 0) are given explicitly by 



a 




%W^-^\c,ayUl^\c,a) 



(2.7) 



The vacuum boundary condition usually^ corresponds to (1, a) = 
(1, 1). The full (r, a) boundary weights are then given by the ac- 
tion of an r-type seam onto the (1, a)-boundary weight. The 
double row seam is given by two regular r-seams sharing the 
same boundary field ^, placed on top of one another, with the 
same spectral parameters as bulk faces appearing in the double 
row transfer matrix defined below in (|2.11| ): 



(r,a) 



.7 
a 



d 6 



b /3 




-tS 



(l,a) 



(2i 



and the left boundary weights are simply equal to the right 
boundary weights. 



^When extra structure is imposed, like in the superconformal case |2l[ , the vacuum 
of the problem can be more complicated. 
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These boundary weights satisfy boundary versions of the equa- 
tions the bulk faces satisfy. The Generalized Boundary Yang- 
Baxter Equation or reflection equation is 




^1 




(2.9) 



^ ' f f 

which is proved using the GYBE (|1.15 ) and the abelian prop- 
erty {^). 



We refer to |jT6|, |T9[ for the boundary crossing equation. 

Let's state here a property that will be of use later on. By 
equation (p.. 7]), one can fill up the triangle appearing in the defi- 
nition ( |2.7| ) of the (1, a) boundary weight with any local face op- 
erators: they will only contribute through a scalar factor, hence, 



9 + 1 , 



2u+(2q-3)\ 



a 




(2.10) 



2.4 Double row transfer matrix 

The double row transfer matrix is given by two rows similar 
to the one appearing in the torus transfer matrix, with spectral 
parameters u for the bottom one and /i — u — (g — 1)A for the top 
one, where is a fixed parameter and q is the vertical fusion 
level. The boundary condition is not cyclic but fixed by the 
boundary weights 
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"JV+l 



fJj-U 



-{<j-i)a) 



-{<j-i)A,Cj 



''JV+2 



0N+2 



"iV+3 



"JV44 



ajv+i ajv+2 <iiv+3 ajV44 




:2.in 



The GYBE dTTsI) implies that double row transfer matrices with 
the same boundary conditions and boundary fields commute: 



(»'L,aL)l(»',s,C)|{»'-R.afl)^ 



(2.12) 



3 Fusion hierachies 

These transfer matrices fulfill a fusion hierarchy of functional 
equations. The details of these equations do depend on the type 
of matrices but their structure is the same. It stems from local 
properties that they all satisfy. Let's choose a horizontal fusion 
level p, and a fixed boundary condition among those available, 
namely toroidal, cylindrical, with or without seams. Call T\{u) 
the corresponding (p, g)-fused transfer matrix at spectral pa- 
rameter u + k\, for — 1 < g < (? — 2, with Tq ^ and Tq defined 
as 

Toi = 0, Tl = f_J ■ (3.1) 
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where is the usual order- bulk term 

fP{ \ — j l-^gi'^)]^ ^ torus 
h ^'^) - j - /i)]^ , for the cylinder. 

(3.2) 

Then the matrices {T'^) -i<q<g-3 fulfill a hierarchy of func- 
tional equations 

n T\ = V,%f^ Tl' + VJl, Tt' (3.3) 

where $q, Vq and Vq, are (scalar) functions that we are go- 
ing to describe, that account for the contributions of the bulk 
faces, the seams and the cylindrical (1, a)-boundary conditions 
respectively. 

The functions Vq and Vq are trivial in the torus case, Vq = 
Vq = 1, and on the cylinder, they are given by 

y ^ Sq_2{2u - fi)S2q+l{2u - fl) 
^ Sq-i{2u- H)s2q{2u- H) 

and 

y ^ Sq{2u- fl)S2q-l{2u- fi) 

The function $q is the product of order-1 terms coming from 
the seams. As we saw in Section an (r, a)-boundary con- 
dition is constructed from the action of an r-seam on a (1, a)- 
boundary condition, and we count separately the type r seams 
coming from the left and right boundaries. If there are K seams, 
the function $q is given by a product of K similar terms: 

K 

% = l[Mrk,Ck,u). (3.6) 

k=l 

The contribution of an (r, s, C)-seam only depends on r and 
(j)q{l, ^, m) = 1. For 2 < r < ^ - 2, 

0* (r, ^, u) = Sq-riu + ^)sq{u + ^) , foT the torus 
<Plir, ^, u) 0g+^_i(r, -/i - ^, m) , for the cylinder. 

(3.7) 
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contribution of 

More generally, we have the following hierarchy of inversion 
identities which can be proved by induction as in |TT| : 

A, Tl T\ = B, f_Jl \[^kI + C, T^+^ T\-' (3.8) 

k=l 

where, in the torus case the functions Aq = Bq = Cq = 1 are 
trivial, and 

Aq{u) = Sq^i{2u - fi)sq+i{2u - /i) , (3.9) 



Bq{u) = s_i(2m - fi)s2q+l{2u - fl) , (3.10) 

Cq{u) = [Sq{2u - fl)]\ (3.11) 

result from the left and right vacuum boundaries in the cylinder 
case. 

If we further define the normalized transfer matrices 

^ rpq—l rpq+l 

then the inversion identity hierarchy can be recast in the form 
of the following universal thermodynamic Bethe ansatz (TBA) 
functional equation 

tit\ = (3-13) 

In deriving the TBA equation we have used the simple proper- 
ties 

Bqiu)BqitL + X) = Bq^l{u + \)Bq+l{u) , (3.14) 

and 

Cq{u)Cq{u + \) ^ ^g^^^^ 



Aq_l{u + \)Aq+l{u) 

Equations ( p.3|) , (|3.8|) and (p.l3|) give a matrix realization of 
the fusion rules ( |1.9| ) and ( p..lO| ). 
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4 Derivation 

Before we proceed to the detailed derivations of ( |0| ) for the in- 
dividual torus and cyhnder cases, let's study the local properties 
which are common to both cases. 



4-1 Local properties 

Firstly, let's look at how the product Tq is decomposed into 
a sum of two terms Tq~^ and Tq"*"^ up to scalar factors. 

Because of the vertical push-through property, we can disre- 
gard the horizontal fusion projectors and apply them later on 
as a wrapping of the equation. 

The product TqTJ is realized as two transfer matrices stacked 
upon each other, the top one being at vertical fusion level 1 and 
the bottom one at fusion level q. Consider an arbitrary column of 
the torus transfer matrix (C-type seams excluded). In fact, after 
a simple manipulation ( [4.13|) , the product of transfer matrices 
in the cylinder case will be built up of similar columns. There 
is a projector P^+i attached to its bottom part, realizing the 
vertical fusion. The Boltzmann weights of this column can be 
written in terms of Temperley-Lieb operators, 

X^,{v) ...X^^{v+iq- 2)A) Xj (v+iq- 1)A) P,^/X,_i(t;+gA) 

(4.1) 

with j an arbitrary label and v the spectral parameter involved 
in that particular column, for example v = u — kX for a typical 
bulk face and v = u + ^ — kX for a face in an r-seam. Because an 
s-type seam is the braid limit of an r-type seam, we don't lose 
any generality in considering only r-type seams. It is easy to 
see that the following arguments can be applied also to C-typs 
seams and that their contribution is trivial. 

We duplicate the projector and insert between its two copies 
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the following identity 
1 



S, 



iS,X,{X) + X,{-qX)) = I 



(4.2) 



v+qX 




p9+l 
J + 1 






^-1)A 



+ 




\ 



(4.3) 



pq + 2 



(4.4) 




The projector -P/^^ in the second term of ( [4.4|) is obtained by 



the definition ( |1.4|) . Thus, this term gives us the column which 



appears in Tq'*'^ in the functional equation (|3.3|). By pushing the 



projector through horizontally in the product of transfer matri- 
ces, we can make it appear in between all columns and because 
of the cyclic boundary condition (and a similar argument in the 
cylinder case), we finally obtain a term which is proportional to 



rjiq+l 

• 



We are now going to prove that the first term of ([4.4| ) yields 
a term proportional to T^^^. 

The product Tq involves a whole row of columns such as 
the LHS of (|4.3|), hence the columns of Boltzmann weights occur 
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with a fusion projector P^+i between each of them. We can use 
the push-through property ( |1.13| ) to remove all but one of these 
projectors one by one: 













A 


v+ 

^9-1) A 


A 




V 

\ 


p<3 + l 


T 






I 


V 


V 


V 




v+ 

^9-1) A 




pq+1 






v+q\ 



(4.5) 



So for cyclic boundary conditions, the projector on the left of 
the +A face in the first term of ( [4.4|) can be discarded and the 
one on its left will be the only remaining projector in the row. 
We will see that the same argument is also valid for cylindrical 
boundary conditions. 

We now make use of the contracting property of the local face 
projector Xj{X): 




SiiU sii-u 




(4.6) 



so that the two faces in the top rows of the first term of 
collapse into a scalar under the propagation of the contractor: 
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(4.7) 
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and the newly appeared contractor further collapses the top two 
faces of the next column on the left. Reapplying the procedure to 
the rest of the columns on the left and using the cyclic boundary 
conditions, we finally collapse all of the top two rows. What is 
left is a scalar contribution —Sq{u)sq-2{u) for each column at 
spectral parameter u, a row of faces with spectral parameter A 
and the local face operator Xj{\) at the right of the projector 
P''+^ But the Boltzmann weight of a face at spectral parameter 
A is simply 

-M^S,,, (4.8) 



W 



A 



hence the top row disappears, leaving 



'H-l "i+-2 




(9-2) A 




p9+l 




(9-2)A 



ai-i ai 



We decompose further the projector P'^^l and sum over di to 



(4.9) 



get the shorter fusion projector P^: 



s 



fd^ 

^6, 




di ^(9-1) 



q+1 



p9 



> ai ■ 



(4.10) 

so that (|4.7| ) reduces to the product of the scalar contribution 
for each column times the matrix valued function Tg~^. 

We now give the details of the contribution of each column 
for each boundary condition. 
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4-2 Functional equation on the torus 

Each horizontally p-fused bulk column in Tq brings a scalar 
factor of s^{u) when collapsed by the contractor. Hence the 
bulk faces contribute to the Tq~^ term as 

/,» = KHr- (4.11) 

The contribution of this same column to the Tq"*"^ term comes 
from the removal of the common scalar factors which appear in 
the process of vertical fusion of the top (p, l)-fused face with 
the larger {p, g)-fused face, yielding a {p,q + l)-fused face. The 
result is fg_i{u). 

Likewise, an r-type seam contributes in the same proportion 
but with a shift in the spectral parameter and an adjustment 
in the common scalar factors, yielding ( p.7| ). It is easily checked 
that the braid limit of such a factor simply vanishes, hence the 
s-type seams don't contribute to the TEA equation and the 
same holds for C-type seams. 

4-3 Functional equations on the cylinder 

As we discussed in Section an (r, a) boundary is the com- 
bination of an r-seam and a (1, a) boundary so we restrict our- 
selves to (1, a)-boundary conditions. 

In the cylinder case, the product Tq of double row transfer 
matrices is realized as four layers of rows and a typical column 
is the stack of two (1, l)-faces on top of two (l,g)-faces, re- 
spectively at spectral parameters fi — u + ^ — qX, u + ^ + qX, 
/i — M + ^ — (g— 1)A and u + C, where ^ = —k\ for a usual 
bulk term involved in a horizontally fused face. Consider the 
following inversion relation ( |1.16D 

xJ^(2m - /X + (2g - l)X)Xf{-2u + /i - (2g - 1)A) = 

(4.12) 
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It follows from the GYBE (gJ^) that 
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2-29 
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q\ 


\ / —2u+fi / 




-l)A\/-(2q-l)A / 




-(g-l)A 








4.13) 



As the rows on the right are of the same structure, the face 
X'^^{—2u + yU — (2g — 1)A) can push through rightward all the 
way to the right boundary. Similarly, its counterpart X^'^{2u — 
yU + (2g — 1) A) can push through leftwards all the way to the left 
boundary. Because of equation ( p.lO|) , these rectangular weights, 
after a crossing symmetry (|1.2CI| ), simply agglomerate into a 
larger boundary (minus the larger projector). For a {l,afj)- 
boundary condition on the right hand side, it reads: 



— TT 

War 



s!_\(/x-2m) 




(4.14) 



In the bulk, we are now in a similar configuration as in the 
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torus case, simply the columns are doubled. Hence the same 
technique applies provided the push through of the contractor 
X{X) and the longer projector P'?+2 behave as expected on the 
boundary. This is what we discuss next. 

It is clear from ([4.14|) that the push-through property of pro- 



jectors is still satisfied on the boundary, the projector P^+i com- 
ing from the bottom q rows of the product can go up through 
the right boundary and get back to the q intermediate rows, 
lower part of the top half. Similarly, coming from the right of 
these rows, it can go down the left boundary to the lower q rows. 
Therefore, the term proportional Tq^^ proceeds in exactly the 
same way as in the case of the torus. 

Consider now the term Tq~^. We need to understand the ac- 
tion of Xj{X) on the (1, aji) right boundary. 

Similarly to ([4.6|) , we have 




Si(m)si(-m) 




di h. 



(4.15) 



So that when the contractor acts on the (1, aji) right boundary, 
we get 





CLR 

(4.16) 
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where we have used the crossing symmetry and then the abelian 
property ( |1.17| ) to interchange the parameters A and 2u+2q\— n 
between two face weights. Then, we apply ( [4 .61 ) to collapse the 
faces weights inside the (l,a/j) boundary 

bjV-2 fej v-1 /)jy- 



i]V-l, 



5L-3(2« - lj)sl~\{2u - fi) 




(4.17) 

where we use the identity ( |1 . 7| ) to eliminate the face weight with 
parameter 2u + 2qX — fi. 

We can see in ( ^T7|) , that the right (1, an) boundary of T^T} 
is contracted into a smaller under the action of the con- 

tractor Xj{X). We then continue to collapse the face weights in 
the bottom half of the transfer matrix with the contractor. 

For the {l,aL) left boundary, the contractor Xj{\) acts from 
the bottom. We rotate the whole diagram by half a turn and 
we apply the same technique, and we get the following scalar 
factors 

(-1)^S2,-3(2m - I2)sl-\{2u - ^i)sl-\{2u - ^i) . (4.18) 

Finally, the contractor can go back to the top row from the 
left, hence the rest of the proof proceeds as previously. 

Collecting the different contributions for all the columns, which 
come in pairs, for the lower and the upper halves of the prod- 
uct Tq T\, with spectral parameters u + ^ and ^ — u — qX + ^ 
respectively, one gets the result listed in Section 
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5 Discussion 



In this paper we have derived the TBA functional equations for 
critical lattice models using simple fusion projectors. We point 
out, however, that the very same functional equations can be 
derived off-criticality by using the methods of [0 . This applies. 



for example, for the A and D models which admit elliptic solu- 
tions to the Yang-Baxter equations. 

We conjecture generally that the form of the TBA functional 
equations are universal for all integrable lattice models associ- 
ated with rational CFTs and their integrable perturbations. In 



particular, we expect the known forms |T0[ of these equations 



to apply to all integrable boundary conditions. 
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